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DIURNAL VARIATION OF THE ION TEMPERATURE* 

K. K. Mahajan** 

Arecibo Ionospheric Observatory , Arecibo , Puerto Rico and 
Goddard Space Flight Center, Greenbelt, Maryland 


ABSTRACT 

Results on the diurnal variation of the ion temperature , obtained 
from the Thomson scatter spectra at Arecibo, in the height range 
250-475 km are presented. At 250 km where the ion temperature 
<Ti ) is expected to be nearly equal to the neutral temperature (T n ), 
T; is generally observed to show diurnal maximum at 16 hr L.T. 
and diurnal minimum at 4 hr L.T. , with no appreciable change in 
these times with solar activity. The average day to night T, ratio 
at 250 km is found to be 1.45. These results are in apparent con- 
flict with the satellite drag analysis in which T n is found to have a 
diurnal maximum at 14 hr L.T. and day to night temperature ratio 
around 1.30. At 475 km where Tj is significantly above T n , the di- 
urnal behavior of T ; is observed to be different on different days. 
The day-to-day and hour-to-hour variations in T ; at this height are 
found to be due to the changes in the electron density and electron 
temperature and seem consistent in terms of an energy budget that 
includes heating of the ion gas by the hotter thermal electrons and 
cooling by the neutral atmosphere. 


♦This work was initiated at the Arecibo Ionospheric Observatory, Arecibo, Puerto Rico 
and completed at Goddard Space Flight Center 
**NAS-NRC Postdoctoral Resident Research Associate 
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I. INTRODUCTION 


It was shown by Hanson (1963) that for altitudes below 300 km the energy 
transferred to the ion gas by the hotter electron gas would be rapidly lost to the 
neutral atmosphere and consequently the ion temperature (T; ) would be very 
nearly equal to the neutral temperature (T n ) . At higher altitudes , due to sub- 
stantial decrease in the neutral gas concentration, the ion gas can not efficiently 
loose the thermal energy (obtained from the electron gas) to the neutral gas without 
the existence of an appreciable difference between T s and T n . Hence a charac- 
teristic transitional behavior of T. , increasing from T at 250 km towards the 
electron gas temperature (T e ) above 700km is expected (Hanson, 1963; Geisler 
and Bowhill, 1965; Dalgarno etal. , 1967; Banks, 1967). As such, experimentally 
observed Tj at altitudes below 300 km can provide information on the diurnal be- 
havior of T n and at higher altitudes can be used to examine the energy budget 
dealing with the heating and cooling of the ion gas. The purpose of this commu- 
nication is to present results on the diurnal variation of T; in this context. The 
ion temperatures to be reported are obtained from the Thomson scatter spectra 
observed in the height range 250 to 475 km at Arecibo during the period October 
1965 through March 1967. 

II. THE METHOD 

The Arecibo facility and its operation for the ionospheric Thomson scatter 
measurements has been discussed by Gordon and Lalonde (1961) and by Carlson 
(1965) . The power spectra of the Thomson scatter signal provide information 
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on the electron and ion temperatures. In the present analysis, these have been 
collected with a 500ps pulse length (resulting in 75 km smearing) by using the 
digital autocorrelation technique developed by Perkins and Wand (1965). Tem- 
peratures from the spectra are derived by two methods. In one of these the half 
power width and the wing to center ratio of the power spectra are used (see e.g. , 
Evans and Loewenthal, 1964). In the other, temperatures are determined from 
the first zero crossing and the depth of the first null of the autocorrelation func- 
tions (Perkins and Wand, 1965). A pure 0 + ion is assumed in interpreting the 
spectra. Temperatures used are the average of the two methods for any one ob- 
servation. At times, in particular during the night, when temperatures by either 
of the two methods were not available (due to low signal to noise ratio or due to 
some inherent difficulties with the method) , these times were not considered in 
the analysis. 

m. THE RESULTS 

Figures 1 and 2 are the diurnal plots of the ion temperature at heights of 
250 and 325 km for several days during the period October 1965 through March 
1967. The observations for each experiment, in general, start around 16 hr. 
L.T. on one day and end at about 20 hr. L.T. next evening and have a time res- 
olution of about 30 minutes. Since Tj at 250 km is expected to be very nearly 
equal to T n (Hanson, 1963; Geisler and Bowhill, 1965; Dalgarno et al. , 1967: 
Banks, 1967), Fig. 1 should basically represent the diurnal behavior of T n . 
While Tj can be more than T n by a lew tens of degrees at 325 km, the difference 


between T, and T n is not expected to change significantly with time during the 
day (Dalgarno et al. , 1967). Nisbet (1967), in fact, has calculated T n at 250 and 
325 km from the Thomson scatter measurements of T e , Tj and electron density 
(N e ) at Arecibo, by using the ion continuity equation. He has found T n to be 
equal to Tj during the night at these heights. During the day he has estimated 
Tj -T n to be of the order of 25°K at 250km and 100°K at 325km with no appre- 
ciable variation in this value with the time of the day. As such Fig. 2 should 
also represent gross features of the diurnal behavior of T n at these heights. 

The day-to-day changes in the value of Tj (at fixed times) are due to changes in 
the solar activity, as demonstrated in an earlier communication (Mahajan, 1967). 
The ’short term' hour-to-hour changes are generally within the uncertainty of 
the measurements (estimated to be about 5%) . The diurnal maximum seems to 
occur most often around 16 hr L.T. (spread between 14 hr and 18 hr on other 
days) and is in disagreement with the 14 hr maximum evolved from the satellite 
drag data (see e.g. Jacchia, 1965). The time of the diurnal minimum is approx- 
imately before local sunrise. 

The times of the occurrence of the diurnal maximum and minimum are more 
clearly demonstrated in Figs. 3 and 4, where the quantity Tj /T min (Tj is the ion 
temperature at any time and T min the minimum nighttime ion temperature for 
that day) is plotted against time for a height of 250 km for several geomagneti- 
cally undisturbed days . Due to minor fluctuations in the nighttime temperatures , 
some smoothing has been done to obtain T min . Figure 3 includes observations 
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on October 5-6, October 27-28, 1965 and January 31- February 1, February 28- 
March 1, 1966 during period of low solar activity. Figure 4 includes observa- 
tions on August 16-17, November 12-13, 1966 and March 3-4, March 23-24, 

1967 during periods of moderate solar activity. The daily and the 27- day average 
(centered at the day of observation) values of the 10.7 cm solar flux, reported by 
the National Research Council, Canada are given in Table 1. It can be noted 
from Figs. 3 and 4 that the average times of the diurnal maximum and minimum 
are respectively 16 hr L.T. and 4 hr L.T. , both during the low and moderate 
solar activity. The average ratio T max /T min is 1.45 ±0.15 and is well above 
the value of 1.28 deduced from the satellite drag data (Jacchia, 1964), but is 
close to the COSPAR model (CIRA, 1965). 

A comment that the 75 km smearing resulting from the long pulse length for 
the Tj measurements at 250km is in order. With the lower edge of the pulse in 
a region of heavier ions (NO + and 0 2 + in small quantities) and the possible in- 
crease of Tj with height in the region occupied by the pulse (i.e. , 212.5 to 287.5 
km), some errors in Tj measurements might be expected. While the errors due 
to temperature increase within the pulse length are expected to be negligible, 

(for nearly constant temperature gradient) , as the measurements correspond to 
the center of the pulse, some underestimation in Tj is possible due to the pres- 
ence of heavier ions near the lower edge of the pulse. It is however interesting 
to note that the Tj measurements are not very sensitive to the presence of heavier 
ions in small quantities. Wand (private communication) nas made some estimates 
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of the error introduced due to the presence of heavier ions by comparing the 
theoretical autocorrelation functions with those experimentally observed from 
the 'double pulse' technique which gives 15 km resolution in height. He has found 
that the presence of 10% of the heavier ions throughout the pulse length under- 
estimates T ; by no more than 2% and by less than 5% with 20% of the heavier 
ions . As the heavier ions for our long pulse measurements are present only in 
lowermost part of the pulse, the errors in this case are expected to be much 
smaller. Ion composition measurements during March 1966 (Brinton, private 
communication) , a period closest to our observations , indicate about 20% of the 
heavier ions at the lower edge of the pulse, decreasing to about 1% at the upper 
edge. From these measurements one would estimate about 5% as an average 
amount of the heavier ions present throughout the pulse length. This would re- 
sult in an underestimate of less than 2% in the T; measurements. It is, how- 
ever, pointed out that these errors are not expected to anyway affect the diurnal 
behavior of T s at 250 km, as the errors will be systematic. Reasonable diurnal 
variation in the ionic composition will probably not significantly alter the results , 
Support to this conclusion also comes from measurements at 325 km where 
the problem of long pulse is less serious (due to pure 0 + region). The diurnal 
behavior of T ( at 325km follows very closely the T ; behavior at 250km. Also, 
as the T n calculated by Nisbet (1967), from the ion continuity equation, at 250km 
is in reasonable agreement with that calculated at 325 km, our confidence in Tj 
measurements at 250 km all the more increases . 
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The diurnal behavior of the ion temperature at a height of 475 km, where T s 
is observed to be significantly above T n , is examined in Fig. 5. Here the day- 
time T ; is plotted against time for three observations. No nighttime data is used 
due to uncertainty in the ionic composition at this height during the night. The 
simultaneously measured T e and N e values for these observations are also plotted 
in Fig. 5 for comparison. The diurnal variation of T, is not quite the same on all 
the three days and is very different from the Tj variation at 250 and 325 km shown 
in Figs. 1 and 2. 

The day-to-day and hour-to-hour changes in T ( can be examined in terms of 
the energy budget dealing with the heating and cooling of the ion gas . The heat 
given to the ions by the hotter thermal electrons is given by: 

7.7 x 10- 6 N e Nj (r e -Ti) 

Qei = d) 

Aj T e 3/2 

where A; is the ionic mass in atomic mass units and Nj is the ion density. As- 
suming oxygen atoms to be the major coolant for the ions at 475 km, then the heat 
given by the ions to the neutral gas is (Banks , 1967): 

Qin =2.1xl0- 15 N e n(0) (T, -T n ) -y/Tj +T n (2) 

n(0) being the atomic oxygen concentration 

If Q ei = Q in and assuming N e = Nj , then from equations (1) and (2): 

Ti-T n / 

T e 3/2 7Ti +T n oc N e 

T e -Ti (3) 
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Thus one should expect a linear relation between N e and the quantity on the 

R-H-S. of equation (3). This is examined in Fig. 6 where the experimentally ob- 
T- -T 

L i x n 

served N e and T e 3/2 «/Tj +T n are plotted against each other for the three 

Te-Tj 

days shown in Fig. 5. Afternoon values between 12-16 hours have been used and 

T n has been taken as the T { measured at 250km. There is a positive correlation 
T- -T 

l n 

between N e and T e 3/ VTj +T n , thus indicating a reasonable agreement be- 

T -T- 
1 e L i 

tween the experimental observations and the theory. Ideally one should expect 

the points to lie on a straight line passing through the origin. This departure 

from the theory could be due to a variety of reasons; the major one being the ac- 

T- -T 

i n 

curacy of the quantity . Although T e and T, are estimated to be correct 

T -T- 

A e A 1 rp 

x i A n 

within ±5%, these limits produce larger error in the ratio Variations 

T e -Tj 

in the number density of atomic oxygen would also add to the scatter. 

IV. DISCUSSION 

From above it is clear that the diurnal behavior of Tj at 250 km (and 325km), 
where T ; is expected to be equal to T n , is in conflict with the neutral tempera- 
ture results derived from the satellite drag analysis. Assuming that the T; 
values presented here are correct, this suggests that either Tj is not equal to 
T n at 250 km or diurnal behavior of T n deduced from satellite drag data is not 
reliable. As already pointed out, Nisbet (1967) has solved the ion continuity 
equation by using Thomson scatter measurements of Tj , T e and N e at Arecibo 
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and has calculated that Tj is very nearly equal to T n at 250km. This would mean 
that the diurnal behavior of T n derived for satellite drag data is probably ques- 
tionable. It is, of course, well known that the satellite drag measurements do 
not provide continuous data throughout the day at one location. (See Priester et 
al. , 1967, for a survey on Satellite drag analysis.) To construct a diurnal pic- 
ture of the temperature, several days data corresponding to various latitudes, 
longitudes , altitudes , solar time , solar activity and magnetic activity might have 
to be used. The diurnal variation of the temperature is then evolved by adjusting 
the derived temperature for variations in the other parameters under suitable 
assumptions . 

At this point it would, probably, be relevant to compare our experimental 
observations with some theoretical results. Harris and Priester (1962) have in- 
vestigated the diurnal variatk of the atmospheric temperature by integrating 
the time dependent heat conduction equation. They found that with EUV as the 
only heat source, the diurnal maximum should occur at 17 hr and have a maxi- 
mum to minimum temperature ratio of 2.2. Although the time of the diurnal 
maximum is close to the time of T, maximum (at 250 and 325 km) shown by our 
measurements, the amplitude is much higher. In the theoretical work of Harris 
and Priester, however, only the vertical expansion of the atmosphere in a column 
was considered. As is now, generally, believed that the horizontal movement of 
the air, giving rise to lateral heat transport, can significantly affect the phase 
and the amplitude of the diurnal maximum (Volland 1966, 1967), any value of the 
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phase and the amplitude can be obtained by choosing suitable phase relations be- 
tween the wind and the temperature. 

V. CONCLUSIONS 

In conclusion it can be said that the diurnal behavior of T; at 250km (and 
325 km) is significantly different from the T n results derived from satellite drag 
analysis. As T ; has been calculated to be nearly equal to T n at 250km by Nisbet 
(1967) , it is concluded that the diurnal behavior of T n evolved from satellite drag 
observations is questionable. The 16 hr L.T. diurnal maximum observed in the 
Tj measurements supports the recently reported Thomson scatter results of 
Carru et al. (1967) and Nisbet (1967). While the diurnal variation of T ; at 
475 km can be different on different days, it seems consistent with the theory 
that assumes heating of the ion gas by the hotter electrons and cooling by the 
neutral atmosphere. 
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Table 1 


Date 

S 10.7 

Daily Values 

® 10.7 

27 Day Average 

Oct. 5-6, 1965 

91, 85 

81, 81 

Oct. 27-28, 1965 

78, 77 

71, 71 

Jan. 31-Feb. 1, 1966 

78, 80 

87, -86 

Feb. 28-Mar. 1, 1966 

86, 81 

81, 81 

Aug. 16-17, 1966 

93, 94 

104, 104 

Nov. 12-13, 1966 

129, 129 

113,.. 114 

Mar. 3-4, 1967 

203, 215 

162, 163 

Mar. 23-24, 1967 

151, 163 

148, 148 
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LEGENDS TO THE FIGURES 


Fig. 1 Diurnal variation of Tj at 250 km for several days. The serial numbers 
are the codes for the experiment. The zero for T, corresponding to 
each experiment is indicated on the vertical scale. 

Fig. 2 Diurnal variation of T, at 32£> km for several days. The serial numbers 
are the codes for the experiment. The zero for Tj corresponding to 
each experiment is indicated on the vertical scale. 

Fig. 3 Diurnal variation of the quantity Tj /T m j n at a height of 250 km for 4 
days during low solar activity. 

Fig. 4 Diurnal variation of the quantity Tj/T m j n at a height of 250 km for 4 
days during moderate solar activity. 

Fig. 5 Diurnal variation of Tj , T e and N e at 475 km for three days. 

T- -T 

i x n 

Fig. 6 A plot of N e against T e 3 '“^/Tj +T n for observations shown 

T -T- 

x e x l 

in Fig. 5. Times between 12-16 hours have been used. 
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